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We report a novel mode of oscillatory crack propagation when a cutting tip is driven through 
a thin brittle polymer film. Experiments show that the amplitude and wavelength of the oscil- 
latory crack paths scale linearly with the width of the cutting tip over a wide range of length 
scales but are independent of the width of the sheet and of the cutting speed. A simple geomet- 
ric model is presented, which provides a simple but thorough interpretation of the oscillations, 
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Despite a long history of research into the field of fracture, many puzzling issues remain unsolved. An interesting 
problem concerns the direction of propagation of the crack tip: when a glass breaks, can the shape of the resulting 
£JQ. pieces be predicted? Recent well controlled experiments have yield a variety of interesting behaviours that are a 
challenge to existing theoretical formulations. An oscillatory instability in dynamic cracks was recently observed in a 
<^ . pre-tensioned thin rubber sheet 6] whose mechanism is still unclear. Another example is the controllable quasi-static 
propagation of oscillatory cracks in a thin strip of glass submitted to a thermal field QI^OI which, despite its apparent 
simplicity, has been stimulating much theoretical study 0, H • 

Here we report novel results on oscillatory fracture paths in a new experimental context: an object, which we 
denote by cutting tip, is perpendicularly driven through a thin polymer sheet held along its lateral boundaries, and 
progressively cuts the material as it advances. For large enough cutting tip widths, the crack follows a well defined 
and highly reproducible oscillatory path that spans a wide range of length scales (Figure Q). In fact, even doing the 
experiment by hand yields surprisingly regular patterns. The experimental observations of oscillatory motion in this 
new geometry present a challenge, from a fundamental point of view, to our understanding of crack propagation in 
thin sheets. Moreover, these ideas should have practical applications since the precise cutting of brittle thin sheets is 
common in industrial manufacturing. 

We perform well controlled experiments with thin sheets of different polymeric materials in a range of thicknesses 
and investigate the dependence of the resulting fracture paths on the width of the cutting tip. We also propose a 
simple model that correctly reproduces the behaviour of the fracture and illustrates the geometrical origin of the 
£h . oscillation. 
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Experimental oscillatory fracture paths 



A schematic diagram of the apparatus is presented in Figure It consists of a thin flat sheet (dimensions ranging 
from 6 x 60 mm to 12 x 50 cm) clamped along its lateral boundaries and mounted on on a linear translation stage. 
This stage was driven at constant speed v towards a fixed object, the cutting tip, which could have either rectangular 
or cylindrical profile with a variety of widths (0.05 mm < w < 60 mm). A camera was mounted directly above the 
apparatus such that the propagating crack was imaged in the cutting tip's frame of reference. 

The sheet was initially prepared with a notch on one of its side boundaries to position and initiate the crack. Both 
bi-oriented polypropylene and cellulose acetate thin sheets were investigated, with thicknesses ranging between 25 and 
130 fim. The sheet's Young's modulus and fracture energy were measured to lie within the ranges E = 1-2 GPa and 
7 = 2-5 kJ/m 2 , respectively. Although they are made of polymers, these materials are brittle as they were severely 
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FIG. 1: Right: edge of the sheet seen from above (polypropylene 27 fim thick) for two diameters of the cylindrical cutting tip: 
w = 31mm and 5mm. (a) Raw experimental snapshots of the cracked film edge and (b) after rescaling according to indenter 
radius; (c) numerical simulation. 
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FIG. 2: A cylindrical cutting tip is forced into a clamped sheet with a notch. The cutting tip is displaced at constant speed 
v with respect to the plate, leading to an oscillatory crack path. Initial configuration of experimental setup (left), typical 
configuration during fracture (center) and geometrical 2D model (right). The strip is flexible in the smallest convex S (light 
grey region) containing the crack path P, where it can easily bend without stretching. Penetration of the cutting tip beyond 
this region induces tensiles stresses in the plate. Our simplified 2D geometrical model takes into acount this important feature 
through the angle a at the tip of the active zone L shown in dark grey: propagation is set to take place in mode III loading 
(inset) along a direction given by /3 when a exceeds a critical value a c . 




(a) (b) (c) (d) (e) (f) 




FIG. 3: Similar sequential processes of the crack propagation over a single period (top: experimental cutting tip in black; 
bottom: corresponding frames in simulation with a = .153 and f3 = 7r/2, crack path in dashed line). Frames (a-c): the left 
hand side of the crack edge is active and the crack tip advances leftwards. In between (a) and (b), the cutting tip crosses the 
upper part of the crack curve and a new contact appears. A kink then takes place, followed by a dynamical propagation (double 
arrow), both in the experiments and in the simulation, up to (d) which is the mirror image of the initial frame (a). The second 
half-period up to (f) is similar. In the simulation, the kink occurs because the active zone (dark grey area) is made up of two 
components, which rest on either side of the indenter. Change of crack direction between (c) and (d) occurs when the angle a 
corresponding to the right hand component of the active zone becomes larger than that of the left hand side. Movies available 
at http: //www. 1mm. jussieu. f r/platef racture 

stretched when processed into thin sheets. As a result, they undergo minimal plastic deformation during fracture 
propagation but, being thin, can sustain large bending without crack initiation. This explains why they are widely 
used in the packaging industry (resistant but easy to tear once a notch is started). The oscillatory paths discussed 
below were not observed in ductile materials. 

As the thin sheet is forced through the fixed tip, the material is cut, leaving behind a well defined and highly 
reproducible oscillatory fracture pattern, some examples of which are shown in Figure ^ As one would expect, for 
thin enough objects the path is straight. The study of the transition from straight to oscillatory paths will be presented 
elsewhere Q and here we focus on the oscillatory regime for large cutting tips, that is well above threshold. 

Although these oscillatory patterns are reminiscent of those observed in thermal quenching experiments 0, U , 
an important difference is that the oscillation mechanism here arises from a coupling with the large out-of-plane 
deflection of the thin sheet, as seen in the experimental frames in Figure 01 
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A simple geometrical model 

We believe that the origin of the oscillations in our experiment lies in the mechanical properties of thin elastic 
sheets, and in their connection with the geometry of surfaces. Under external constraints, thin elastic sheets bend 
in order to avoid stretching, since bending energy is small compared to stretching energy 7]. In order to conserve 
their in-plane lengths and remain unstretched, thin sheets tend to adopt developable shapes, that are the reunion 
of straight lines, generatrices. Because the lateral edges of the plate are held, these generatrices may not cross the 
clamped boundaries of the plate and can only end up on the crack path. We thus define a "soft" zone in the plane 
of the elastic plate by considering the reunion of all the segments ending on the crack path. This is called the convex 
hull S of the crack path (see Fig. [2J right). In 5, the sheet can deform by pure bending and can therefore easily 
accomodate the presence of the cutting tip. However, when the cutting tip penetrates beyond 5, in a region which 
we call "hard" , the sheet undergoes stretching which eventually leads to crack propagation. 

In order to determine the direction of propagation in our thin sheets we have performed a complementary ripping 
(mode experiment by systematically curving and pulling one side of an initial notch. Crack propagation was 
observed to always follow a well defined and highly reproducible angle (3 with respect to the generatrix of the fold. 
This angle (3 depends on the material and, since these materials are anisotropic, on the orientation of the ripping, but 
we consistently measured an angle close to tt/2 to within 20%. So, to first approximation, we take (3 = tt/2. 

The propagation criterion in our model is the following. We consider the "large" hull A, defined as the convex hull 
of the crack tip plus the horizontal section of cutting object (black disk in the Fig. 2c) altogether. The so called active 
area L (in dark grey) computed as the set difference L = A \ S of the large and the soft hull is the region of the hard 
zone where penetration of the cutting tip induces tensile stresses. When this "active" area includes a cone with angle 
a > a c at the crack tip greater than a critical value (see Fig.Efc), the crack is set to propagate instantaneously in the 
direction given by the angle /?, until a decreases back to a c (note that a indeed decreases during crack advance since 
the soft zone expands). The angle a c results from a balance between in-plane stresses in the plate and the energy 
needed to advance the crack: it can in fact be computed from the mechanical and fracture properties of the film. We 
refer to Ref. for a more detailed presentation of this geometrical fracture model and its rational derivation from 
elementary principles of fracture and plate mechanics. 

Experimental and numerical results 

A numerical simulation of this geometrical model yields oscillatory crack paths over a wide range of values of the 
parameters a c and (3. In Fig. [I] we compare typical computed (right) and experimental (left) crack paths. The 
non-sinusoidal crack patterns are in remarkably good agreement. Their wavelength and amplitude depend weakly on 
the parameters a c and (3 of the model and their overall shape is almost independent of them. 

The sequential processes of the crack propagation over a single period is presented and discussed in Figure 01 for 
both the experiments and numerics. The numerical simulations were run with parameters a = .153 and (3 = tt/2. 
The evolution of the crack paths demonstrates strong similarities between the experiments and the simulation based 
on our model. 

Our geometrical model includes a single length scale, the size of the cutting tip, and does not take into account its 
driving speed. In fact, we have performed experiments for 0.06 mm/s < v < 64 mm/s and found no influence of speed 
of the patterns. We have also checked experimentally, for a fixed tip size w = 2 mm, the pattern to be independent of 
the sheet's width, D, within the range 3.3w < D < 21w, in agreement with the predictions of the model. This feature 
is a priori surprising since one could have thought of D as one of the natural length scales. Having determined that 
there is only one length scale w in the problem, one expects both the wavelength and amplitude of the oscillations to 
scale linearly with it. Indeed, as shown on Fig. we experimentally find a linear dependence of the amplitude and 
wavelength on the tip width over almost three decades, for a variety of polymeric materials and thicknesses. This 
supports a purely geometric fracture mechanism in this regime. Note that this linear behaviour breaks down in the 
vicinity of the transition region between straight and oscillatory paths: there, the sheet can no longer be considered 
thin relative to the cutting tip and the model's assumptions are not valid anymore (the plate bending and stretching 
energies become of the same order of magnitude). 

To conclude, we have reported a novel mode of oscillatory crack propagation as a cutting tip of moderately large 
width is driven through a thin polymer sheet. We proposed a geometric model for these cracks, which retains the 
salient mechanical properties of thin elastic sheets: our scenario is based on the coupling between the direction of 
crack propagation and large deflections of the plate. Although very simple, such a model reproduces surprisingly many 
experimental features of the cracks, including their shape and dynamics. The crack path amplitude and wavelength 
were found independent of the driving speed, of the thickness and lateral width of the sheet, and scale linearly with the 
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FIG. 4: Log-log plots of the peak-peak amplitude and wavelength of the oscillations as a function of tip width w for various 
materials and thicknesses (polypropylene 25 to 53 /im: open symbols, cellulose acetate 100 and 130 /im : closed symbols), and 
for both cylinders (om) and rectangular blades (□) : amplitude and wavelength organize on slope 1 lines, on a range of almost 
3 decades. The straight lines represent the best linear fits with slope 1. 



size of the object, in agreement with the model. This explains why this phenomenon is highly robust and insensitive 
to irregularities in the speed and boundary conditions, and can even be observed if the cutting is done by hand. 

The authors wish to thank Mokhtar Adda Bedia for stimulating discussion, and Charles Baroud for putting in 
contact the Paris and Manchester groups who started working on this subject independently. 
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